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We perform an amplitude analysis of B — > ^(1020) K* (892) decay with a sample of about 
384 million BB pairs recorded with the BABAR, detector. Overall, twelve parameters are measured, 
including the fractions of longitudinal /l and parity-odd transverse f± amplitudes, branching frac- 
tion, strong phases, and six parameters sensitive to CP-violation. We use the dependence on the 
Kir invariant mass of the interference between the J p = 1~ and + Kir components to resolve 
the discrete ambiguity in the determination of the strong and weak phases. Our measurements of 
f L = 0.49 ± 0.05 ± 0.03, f± = 0.21 ± 0.05 ± 0.02, and the strong phases point to the presence of a 
substantial helicity-plus amplitude from a presently unknown source. 

PACS numbers: 13.25.Hw, 13.88.+e, 11.30.Er 



The polarization anomaly in vector-vector charmless 
hadronic B-meson decays 0, 0, |, H i, H, Q motivates 
a revision in our understanding of the effective flavor- 
changing b — > s quark transition in B-meson decays. Ex- 
planations of this anomaly led to development of models 
either with physics beyond the standard model 0], new 
weak dynamics [i} , or strong dynamics [l(| ■ 

A vector- vector B-meson decay, such as B — > ipK*, 
is characterized by three complex helicity amplitudes 
Ai\ which correspond to helicity states A = — 1,0, +1 
of the vector mesons. The Aio amplitude is expected 
to dominate ll| due to the (V — A) nature of the weak 



interactions and helicity conservation in the strong in- 
teractions. Experimental results suggest that Ai + i and 
comprise about 50% of the total decay amplitude 
m B -> ipK* [J 0. Recently, the BABAR experiment 
extended the study of the B° —f <pK* Q decay to resolve 
the discrete ambiguity between the Ai + i and j4i_i am- 
plitudes Q. 

We now investigate the polarization puzzle with a full 
amplitude analysis of the B ± — > ipK* (892) ± decay. In 
this paper, we report twelve independent parameters for 
the three B + and three B~ decay amplitudes, six of 
which are presented for the first time. Moreover, we 
use the dependence on the Kir invariant mass of the in- 
terference between the J p = 1~ and + (Kir)^ compo- 
nents [H,[l2,[l3| to resolve the discrete ambiguity between 
the .Ai+i and A%-i helicity amplitudes. 

We use a sample of 383.6 ± 4.2 million T(4S I ) -> BB 
events collected with the BABAR detector [14j at the PEP- 
II e + e~ asymmetric-energy storage rings. The e + e~ 
center-of-mass energy y/s is equal to 10.58 GeV. Mo- 
menta of charged particles are measured in a tracking 
system consisting of a silicon vertex tracker with five 
double-sided layers and a 40-layer drift chamber, both 
within the 1.5-T magnetic field of a solenoid. Identifi- 
cation of charged particles is provided by measurements 
of the energy loss in the tracking devices and by a ring- 
imaging Cherenkov detector. Photons are detected by a 
CsI(Tl) electromagnetic calorimeter. 

The B ± -> ^(1020)isT* ± -> {K+K-^Kn^ candi- 
dates are analyzed with two (Kir)'* 1 final states, Kgir^ 
and A" ± 7T°. The neutral pseudoscalar mesons are recon- 
structed in the final states K s — > 7r+7r~ and ir° — ► 77. 
We define the helicity angle 9i as the angle between the 



direction of the K or K + meson from A* — > Kir (9\) or 
ip — ► K + K~ (O2) and the direction opposite the B in the 
K* or cp rest frame, and $ as the angle between the de- 
cay planes of the two systems 0. The differential decay 
width has four complex amplitudes Aj\ which describe 
two spin states of the Kir system ( J = 1 or 0) and the 
three helicity states of the J = 1 state (A = or ±1): 



d 3 r 



dHidH 2 d<P 



'A JX Yf(Hi,*)Y{*(-m,0) 



(1) 



where Hi — cos 0i and Yj are the spherical harmonics 
with J = 1 for A* (892) and J = for (Att)^ We repa- 
rameterize the amplitudes as Ai±i = (Ami ± Ai±)/y/2. 

We identify B meson candidates using two kinematic 
variables: m ES = [(s/2 + p r • p B ) 2 /A|- — p^] 1 / 2 and 
AA = (Ex Eg - p r ■ p B - s/2)/ A /s, where (A B ,p B ) is 
the four-momentum of the B candidate, and (Ey, Pr) 
is the e + e" initial state four-momentum, both in the 
laboratory frame. We require tues > 5.25 GeV and 
|AA| < 0.1 GeV. The requirements on the invariant 
masses are 0.75 < m K7T < 1.05 GeV, 0.99 < m KR < 1.05 
GeV, Im^ - m K o\ < 12 MeV, and 120 < m 77 < 150 
MeV for the K* ± , ip, Ag, and 7r°, respectively. For the 
Kg candidates, we also require the cosine of the angle 
between the flight direction from the interaction point 
and momentum direction to be greater than 0.995 and 
the measured proper decay time greater than five times 
its uncertainty. 

To reject the dominant e + e~ — > quark-antiquark back- 
ground, we use the angle 9t between the £?-candidate 
thrust axis and that of the rest of the event, and a Fisher 
discriminant T [lij]. Both variables are calculated in 
the center-of-mass frame. The discriminant combines 
the polar angles of the B-momentum vector and the B- 
candidate thrust axis with respect to the beam axis, and 
two moments of the energy flow around the B-candidate 
thrust axis [TBI ]. 

To reduce combinatorial background with low- 
momentum 7T° candidates, we require 7i\ < 0.6. When 
more than one candidate is reconstructed, which happens 
in 7% of events with A^ and 17% with ir°, we select the 
one whose x 2 of the charged-track vertex fit combined 
with x 2 °f the invariant mass consistency of the A^ or 
ir° candidate, is the lowest. We define the 6-quark fla- 
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FIG. 1: Projections onto the variables (a) rriKn, (b) m K ^, 
(c) AE, and (d) 771es for the signal B — ► ^{Ktt)^ 1 candi- 
dates with a requirement discussed in the text. The solid 
(dashed) lines show the signal-plus-background (background) 
PDF projections. 



vor sign Q to be opposite to the charge of the B meson 
candidate. 

We use an unbinned, extended maximum-likelihood 
fit [HIH to extract the event yields and the parameters 
of the probability density function (PDF) "Pj\ The in- 
dex j represents three event categories used in our data 
model: the signal B ± — ► (p(A"7r) ± (j = 1), a possible 
background from B ± — > f (980)K* ± (j = 2), and com- 
binatorial background (j = 3). The superscript k cor- 
responds to the value of Q = ± and allows for a CP- 
violating difference between the B + and B~ decay am- 
plitudes (A and A). In the signal category, the yield 
and asymmetry of the B ± — > </?A*(892) ± mode, n s i g 
and Acp, and those of the B ± — > ^(Att)^ mode are 
parameterized by applying the fraction of <pA*(892) ± 
yield, ^t fe , to n*. Hence, n S ig = raj x /i+ + nj~ x jii - " 



-4cp = ( rl i x A* — n i x A* )/ n sig, and the ^(Att)^ 
yield is nf X (1 — ^t + ) + x (1 — /i~). 

The likelihood A for each candidate z is defined as 
& = J2j,k njVji^i; fx k , C, I), where the PDF is formed 
based on the following set of observables = {Hi, Hi, 
$, rriK-K, m KKi m ES, J 7 , Q} and the dependence 

on fi k and polarization parameters C is relevant only for 
the signal PDF V\. The remaining PDF parameters £ 
are left free to vary in the fit for the combinatorial back- 
ground and are fixed to the values extracted from Monte 
Carlo (MC) simulation 17fl and calibration B — > Dir de- 
cays for event categories j — I and 2. 

The helicity part of the signal PDF is the ideal angu- 
lar distribution from Eq. (fT]), multiplied by an empirical 
acceptance function Q{H\,H2,^) = Gi(Hi) x ^2(^2)- 
Here, the amplitudes Aj\ are expressed in terms of 
the polarization parameters £ = {fa, f±, <pu, <f>±, 80, 
A° CP , A^ P , A0,|, A0_l, AS Q } defined in Table U CP- 
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FIG. 2: Projections onto the variables (a) Tii, (b) H2, (c) 
Q $, and (d) the differences between the Q "J? projections for 
events with Hi H2 > and with Hi H2 < for the signal 
B — > (^(Avr)* candidates following the solid (dashed) line 
definitions in Fig. [T] The step in the Hi PDF distributions 
is due to the selection requirement Hi < 0.6 in the B — * 
ip(K ± TT ) channel. 



violating differences are incorporated via the replace- 
ments in Eq. (TT]) for B + decays: /i^/jx {1+A° CP x Q), 
fa -> /a. x (1 + x Q), 0|| -» (0|| + A0[, x Q), ± -> 
(0_L + tt/2 + (A0_l + tt/2) x Q), and <5 -► (£ + A<S x Q). 

A relativistic spin-J Breit-Wigner amplitude param- 
eterization is used for the resonance masses @, Eg] , and 
the (Ktt)^ nriK-K amplitude is parameterized with the 
LASS function 0. The latter includes the A ( J(1430) ± 
resonance together with a nonresonant component. The 
interference between the J = and 1 (A^7r) ± contribu- 
tions is modeled with the three terms 2TZe{Ai\A^ ) in 
Eq. ([T]) with the four-dimensional angular and ttlk-k pa- 
rameterization and with dependence on /z fc and £. 

The signal PDF for a given candidate i is a joint PDF 
for the helicity angles and resonance mass as discussed 
above, and the product of the PDFs for each of the re- 
maining variables. The combinatorial background PDF 
is the product of the PDFs for independent variables 
and is found to describe well both the dominant quark- 
antiquark background and the background from random 
combinations of B tracks. The signal and background 
PDFs are illustrated in Figs. [T] and O For illustration, 
the signal fraction is enhanced with a requirement on the 
signal-to-background probability ratio, calculated with 
the plotted variable excluded, that is at least 50% effi- 
cient for signal S ± — > < y 9(A7r) ± events. We use a sum of 
Gaussian functions for the parameterization of the sig- 
nal PDFs for AE, toes, and T . For the combinatorial 
background, we use polynomials, except for m-ES and T 
distributions which are parameterized by an empirical 
phase-space function and by Gaussian functions, respec- 
tively. Resonance production occurs in the background 
and is taken into account in the PDF. 
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TABLE I: Summary of results for the B — > (pK* (892)* decay. The twelve primary results are presented for the two decay 
subchannels along with the combined results, where the branching fraction B is computed using the number of signal events 
n s ig and the total selection efficiency e, which includes the daughter branching fractions 0] and the reconstruction efficiency 
£rcco obtained from MC simulation. The definition of the six CP-violating parameters allows for differences between the B + 
and B~ decay amplitudes A and A with superscript Q — — and ±, respectively. The systematic uncertainties are quoted last 
and are not included for the intermediate primary results in each subchannel. The dominant fit correlation coefficients (C) are 
presented, where we show correlations of So with <f>u/<f>± and of A<5o with A<j!>||/A<^x- 



parameter 



definition 



K*(892) ± -> ifg 7 r ± if*(892) ± -> K+iv 



combined 



B 

h 
fx 
<t>\\ - 1" 

<f>± - 7T 
So - 7T 



r/r t 



\A W \ 2 /Z\A 1X \ 



|j4ix|7E|.Aia 
arg^xy/Aio) - 
aig(A 1± /A 10 ) - 
a,rg(A 00 /A 10 ) - 



2 
2 

■ 7T 

- 7T 

- 7T 



(10.5 ± 1.4) x 10~ 6 

0.51 ±0.07 

22+ - 07 
-0 7^+ - 28 
-0.15 ±0.24 
-0.25 ±0.24 



(11.6 ± 1.5) x 10~ 6 

4fi+ 010 

21 + 09 
-0.77 ±0.35 

-0.89«;« 
±0.11 ±0.31 



(11.2 ±1.0 ±0.9) x 10" 6 
0.49 ±0.05 ±0.03 
0.21 ±0.05 ±0.02 
-0.67 ±0.20 ±0.07 
-0.45 ±0.20 ±0.03 
-0.07 ±0.18 ±0.06 



}-58% 

} ±56% 
±37%/ ± 36% 



Acp 


(r+ 


-r-)/(r+ + r-) 


-0.09 ±0.13 


±0.07 ±0.13 


0.00 ±0.09 ±0.04 


A CP 


(ft 




±0.24 ±0.15 


±0.09 ±0.20 


±0.17 ±0.11 ±0.02 


A CP 


(ft 


-a) /(ft + fi) 


±0.12 ±0.31 


+0 41 + - 54 


±0.22 ±0.24 ±0.08 


A0|| 




(0+ - 0jT)/2 


±0.02 ±0.28 


±0.22 ±0.35 


±0.07 ±0.20 ±0.05 


A^x 


(< 


it - 4>1 - tt)/2 


±0.18 ±0.24 


±0.48«;« 


±0.19 ±0.20 ±0.07 


A5 




(tf-*o~)/2 


±0.13 ±0.24 


±0.34 ±0.31 


±0.20 ±0.18 ±0.03 



}-50% 

} ±57% 
±37%/ ± 37% 



72 s jg 
£ 

^rcco 



102 ± 13 ± 6 
(2.53 ±0.13) % 
(22.3 ± 1.2) % 



1171x6 ±7 
(2.59 ±0.17) % 
(16.0 ± 1.0) % 



We observe a nonzero B ± -> ( y 5if*(892) ± yield with 
significance, including systematic uncertainties, of more 
than IOct. The significance is defined as the square root 
of the change in 2 In £ when the yield is constrained to 
zero in the likelihood C. In Table U results of the fit are 
presented, where the combined results are obtained from 
the simultaneous fit to the two decay subchannels. 

We repeat the fit by varying the fixed parameters in £ 
within their uncertainties and obtain the associated sys- 
tematic uncertainties. We allow for a flavor-dependent 
acceptance function and reconstruction efficiency in the 
study of asymmetries. The biases from the finite resolu- 
tion of the angle measurements, the dilution due to the 
presence of fake combinations, or other imperfections in 
the signal PDF model are estimated with MC simulation. 

The nonresonant K + K~ contribution under the ip is 
accounted for with the B° — > foK*° category. Its yield 
is consistent with zero. The PDF shape in this 

category is varied from the resonant to phase-space and 
the yield is varied from the observed value to the extrap- 
olation from the neutral £?-decay mode [B| to estimate 
the systematic uncertainties. Additional systematic un- 



certainty originates from other potential B backgrounds, 
which we estimate can contribute at most a few events 
to the signal component. The systematic uncertainties in 
efficiencies are dominated by those in particle identifica- 
tion, track finding, and Kg and 7r° selection. Other sys- 
tematic effects arise from event-selection criteria, ip and 
K*° branching fractions, and the number of B mesons. 

The yield of the (p(Kir)^ contribution is 571^ events 
with a statistical significance of 7.9a, combining the 
|^4oo| 2 term and the interference terms 2He(Ai\AQ a ), 
which confirms the significant S-w&ve Ktt contribu- 
tion observed in the neutral B-decay mode The 
dependence of the interference on the Kn invariant 
mass [lH, [H| allows us to reject the other solution 
near (2ir — (f>ii,iT — (f>x) relative to that in Table [I] with 
significance of 6. 3c, including systematic uncertainties. 

The (V — A) structure of the weak interactions, helic- 
ity conservation in strong interactions, and the s-quark 
spin flip suppression in the penguin decay diagram sug- 
gest \A W \ > \Ai+i\ > |j4i_i| [11]. This expectation 
disagrees with our observed value of /j> We obtain the 
solution </>|| ~ <f>±_ without discrete ambiguities, which is 
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consistent with the approximate decay amplitude hierar- 
chy \A 10 \ ~ \A 1+1 \ > L4i_i|. 

We find that 0j_ and ^>ii deviate from either 7r or zero by 
more than 3. fa and 2Aa, respectively, including system- 
atic uncertainties. This indicates the presence of final- 
state interactions not accounted for in naive factoriza- 
tion. Our measurements of the six CP- violating param- 
eters are consistent with zero and exclude a significant 
part of the physical region. We find no evidence of CP 
violation in this decay. 

In summary, we have performed a full amplitude anal- 
ysis and searched for CP-violation in the angular distri- 
bution of the i? ± — > {pK*^ decay. Our results are sum- 
marized in Tableland supersede our prior measurements 
in Ref. [l|. These results find substantial A% + i amplitude 
in the B ± — ► ipK*^ decay and point to physics outside 
the standard model or new dynamics [1, \Q, E3] • 
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